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On the Solution of Certain Types of Linear Differential 
Equations in Infinitely Many Variables. 

By Webster G. Simon. 



The main purpose of this paper is to prove the existence of certain types 
of solutions of particular kinds of linear differential equations with periodic 
coefficients in infinitely many variables. As a means to this end the existence 
of exponential solutions is established for certain types of linear differential 
equations in infinitely many variables with constant coefficients. 

The starting point is the existence theorem given by von Koch,* and a 
generalization of Poincare's theorem -f- concerning the development of the solu- 
tions of the differential equations as power series in a parameter fi when the 
functions appearing in the differential equations are themselves power series 
in fi. Then our work is very similar to that of the finite case, the finite 
determinants becoming infinite determinants which together with all their 
first minors converge absolutely. 

The type of determinant used in this paper is more general than the 
normal determinant, and is the following.§ The determinant 

1 + flu, a 12 , a 13 , 

a 21> 1 ~t~ #225 #23> 

ffl 315 fl 32> 1 "T tt 33> 

> > y 

is such that there exist two sets of positive constants 8 1} 8 2 , ■ • ■ , T 1} T 2 , • ■ • , 

00 

which are of such a nature that | «»,- | < SiTj, and S SiTi converges. This 

is a type of infinite determinant given by von Koch,§ which together with all 
its minors converges absolutely. 

Using the methods thus indicated, we find that many of the phenomena of 
the finite systems are carried over into the infinite systems of differential 
equations. 1 1 



* Von Koch, Ofversigt af Kongliga Vetenskaps Akademiens FSrhandlingar, Vol. 56 
(1899), pp. 395-411. 

t Poincarg, " Lea 'MSthodes Nouvelles de la MScanique Celeste," Vol. I, chapter II. 

§ Von Koch, Acta Mathematiea, Vol. 24 ( 1901 ) , pp. 89-122. Hereafter this paper 
will be denoted by K. 

1 1 See Moulton and MaeMillan, American Journal or Mathematics, Vol. 33 
(1911), pp. 63-96. 
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28 Simon: On the Solution of Certain Types of Linear 

I 

The differential equations which we shall consider are of the general type 

00 

(1) tf. = %6ij(t)xi, i = l, 2, • ■ ■ oo, 

7=1 

where x'. is the derivative of x% with respect to the independent variable t. 

Concerning this system of differential equations we shall establish two 
existence theorems. 

Theorem I. Suppose the system (1) satisfies the following hypotheses: 
(ELJ The 6ij (t) are analytic functions of t when 1 1 \ < B. 
(H 2 ) Positive constants, S u S 2 , • • • , T ± , T 2 , • ■ • , exist such that \ 6u (t) \ 
< SiT h \ t | < B, and S 1 T 1 + S 2 T 2 + • • • converges. 
(H 8 ) The Xi{Q) = fa, where the fa are constants such that faT 1 + faT 2 ■ • ■ 
converges. 

Then there exists a unique system of functions 

00 

(2) Xi = Pi + 2(iyP, t — 1,2, - • - oo, 

i=i 

which satisfy the system (1), and which converge for \ t | < B. 

Although this theorem was established by von Koch,* it will be proved 
briefly here because the notation and the results are essential for the later 
parts of this paper. 

From (.Hi) we have 

(3) $ ti (t) = a<°> + a™t + • • • , i, j = 1, 2, • • • oo. 

Upon substituting the series (2) in (1) and equating coefficients, we get 

00 

<3 (1 > = 5 a m fa°\ 



(4) 



00 



00 \ 

(A + l)j8 < ^+ 1) = 22 a%-»p<*», where 0< O) = fa. 

1=1 fc=o ; 

We see that the formal solution is unique. 
* Loo. oit. 
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To prove the convergence of (2), consider the system of differential 
equations 

00 

(5) e — BaT,$ h i=l,2, • • • oo. 

i=l 

00 

Upon substituting & = fli -f 2 yVP in (5) and equating coefficients, we get 

i-i 

oo 



(6) 



3=1 
00 



2 y ™=S i ST j yf\ 



(A + l) y ^+ 1 > = Sa T n <*>, where y< ( 0) — pt. 
j=i 

On comparison of (6) and (4) we see that yW > | j8 (,) |, i,;' = 1, 2, • ■ • oo. 
Hence & dominates Xi for every i. 
From (5) we have 

v ' Si 1 S 2 2 c 

On taking £(0) =0, we have 
(8) £i — pi=>Si£, » — 1,2, •••oo. 

Therefore each equation (5) reduces to 

(9) ?-«+*, 

where — S.T, + fl 2 !F 2 + • • • , K = &T X + &!T 8 + • • • . 

Prom the theory of a finite system of differential equations we know that (9) 

has a unique analytic solution for 1 1 | < R. On combination of this fact 

with (8), we see that the solution (2) of (1) converges when \t\ < R, 

Now we define as a fundamental set of solutions of (1), a set such that 

every solution of (1) can be expressed as linear homogeneous functions with 

constant coefficients of the elements of the set. Then if we denote by $ij the 

elements of the fundamental set, it follows that the determinant of the <j> i} 

converges and is not zero for all J 1 1 < R, and conversely. Furthermore it 

(2$ u (t)dt 
has the form A = A e-'«= 1 ' for all | t | < R, where A is the value of 

the determinant A when t = t . If, for example, we take the solutions defined 

by <t>u(0) =1 and <j>u(0) = °> i¥=j, we see from (8) and (9) that 

] <t>« I ^1 +~± I (««-l) |,and | 4> i} | S ^p| (e«-l) \,i^j. 
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We see that the determinant of the <j>n is absolutely convergent. 

Theorem II. Suppose that the system (1) satisfies the following hypo- 
theses: 

(H x ) The 6ij are expansible as power series in a parameter p., which converge 

for all real \ t | < B if | p, | < p. 

(H 2 ) For p= 0, the $n = an, where the an are constants. 

(H 3 ) The Bn satisfy all the hypotheses of Theorem I uniformly with respect 

to p. if | p, | < p. 

Then the solutions of (1) can be expressed as power series in p, which 
converge for | p, | ^ p < p and for all real | 1 1 < B. 

00 

In Theorem I we showed that Xi — p t = 2 P^P, \ t | < B, \ p, | < p, 

i=i 
where nB M = P ( ">(<8 (0 > • • • , fl'"- 1 '). 

00 

In particular /J' 1 ' = 2 a^W '- The a^°> are power series in p. Hence 

00 00 00 

# 0) <C = V|>*, jSJ" = 2 2#>* which converge for | p. |< p. Then by 
a well-known theorem in the theory of double series, we infer that p (1) = 

00 

2 P^f^, | /* | < p- In a similar fashion we proceed step by step, and show 

that the p { ^ are also power series in p., which converge for | p, | < p. So we 
have 

00 

(10) z« = 2<Mp), t = l,2, • • • oo, 

k=l 

which is uniformly convergent for | p, | fg p < p, and where <#>« (p) = 
^ifc' "I" ^« /* + ' ' ' • Then t>y a well-known theorem of Weierstrass, we 
know that we can rearrange the series (10), and write it in the form 

00 

z 4 = 2 s</> (*)/*>, |p[^p <p, \t\ <B, »— 1,8, • • • oo. 

i-x 

II 

Now let us assume that the $n of Theorem I are constants ; and let these 
constants be denoted by »»,-. We inquire whether the system (1) has a solu- 
tion of the form 

(11) Xi = de at , 

00 

where by (H 3 ) of Theorem I %aijCj converges. Necessary and sufficient con- 
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(12) 



(a u — a) Cj + ia^c., + a ls c s + • • • =0, 
&2i c i + («22 — <*)c 2 + a 2s c s + ' ' ' = °> 



In order to see that (12) can have a solution for which the c% are not all zero, 
we make the transformation a = — A" 1 , assuming of course that a is not zero. 
Then the equations (12) become 



(13) 



(1 -4- A^Jc,. + Aa^c-s -4- • • • «= 0, 
Aa 21 C! + (1 + A« 22 )c 2 + • • • =0, 



In order that the equations (13) have a solution in which the d are not all 
zero, it is necessary and sufficient that A be a root of 



(12) A(A) 



1 + \a t 



Aa, 



A« 21 , 1 -f- Aa 22 , 



= 0, 



Equation (14) is called the characteristic equation. Since the an satisfy the 
hypotheses of Theorem I, it follows that A (A) is an integral function * of A. 
If A is a root of (14), then there exists a finite number r such that 
A(A ) together with its minors of order 1, 2, • • • , r — 1, vanishes; but there 
is at least one minor of order r which is different from zero.f Let the minor 
which is obtained by replacing the elements in the tjth row and the fc x th 
column by unity and the remaining elements in that row and column by zero, 
and so on up to the i r th row and the & r th column be denoted by 



/i 1} i 2 , ■ • ■ , »,s'W 



and suppose that this minor is one which is different from zero. Then the 
system (13) admits the solution I 



(15) c ft = 



(h, • • • *r \ < Xo> , (h,i 2 ■ ■ ■ ir\ «•» , . . , (h, ■ ■ ■ ir\^ 
V^, • • • hr) *h~ T \lCulc ■ • ■ TCr) S," 1 " "^ U 15 ■ ■ • k) S 



V»U ' ' ' KrJ 



<\o> 



depending on r independent parameters fa, • • • & , and admits no others. 
To sum up our results we have the following 



* See K, p. 104. 



t See K, p. 108. 



t See K, p. 109. 
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Theorem III: If the Oij of (1) are constants and if A is a root of 
A (A) = such that all minors of order r — 1, hut not all of order r, vanish 
for A = A , then there exist r independent solutions of (1) of the form (11). 

Now the characteristic equation, since it is an integral transcendental 
function, may have no finite roots, it may have a finite number of finite roots, 
or a denumerably infinite number of roots. In the first case the d are all 
zero, in the second case, a finite number of solutions of the form (11) exists 
and in the last case an infinite number of the form (11) exists. This is the 
case we shall study, for this infinite set constitutes a fundamental set of solu- 
tions as the following discussion shows. 

For each A that is a root of A (A) = there is a solution of (1) of the 
form 

(16) Xij — cije a j t , i, y=l, 2, • • • oo. 

We know that the system (1) admits the fundamental set of solutions <f>ij, 
i,j=l, 3, • • • oo, where <f>u(0) = 1, <£ ii (0)=0, i^j. Consequently 



x h = 5 A^^nc, i = 1, 2, 

7c=l 



00. 



First let us assume that the a,j are all distinct. Then it follows immediately 
that the system (16) constitutes a fundamental set. For consider x^, <j>n, 
i = 1, 2, ■ • • oo, ; = 2, 3, • • • oo. 






=4< i > 



011> 

$215 


0125 ' ' 
0225 ' ' 


5 


y 



Choose the initial conditions such that A w = 1, which is no actual restriction 
since we are merely determining the arbitrary constant which enters, and we 
can multiply our solution by a constant, not zero, when we are through if we 
choose. Thus x iu 4> ih i = l,2, • ■ ■ oo, /=2, 3, • • • oo, constitute a fun- 
damental set. Similarly we show that x il} x i2 , <j> ih constitute a fundamental 
set, where i = 3, 4, • • • oo. Continuing in this manner and passing to the 
limit we see that the x iS , i,j=l, 2, • • • oo, constitute a fundamental set. 
A result of this determination of the arbitrary constant is that c« = 1. 

Next let a 2 = o^ and ay ^ Oj, / = 3, 4, • • • oo . There are two cases 
that may arise here, viz., all of the first minors of A (A) vanish when a=a iy 
or not all vanish. Consider the former case. Then (15) assures us that we 
can solve for the c i} in terms of two of them. On choosing the notation so 
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that (-. 2 J is not zero, we get two independent solutions of the form 

osii — Ci 1 e a i t , xi 2 = c« 2 e a i*, i == 1, 2, • • ■ <x>, 
which with fa}, ;' = 3, 4, • • • oo, constitute a fundamental set; for 



""ll> *"l2J <Pl85 
X 2V a '22> < f , 23> 



A^fa, + a«</> 22 , A^fa t + 1;% tf> 23 , • 



and the A a \ A i2) may be taken to he zero, and A a) = A m = 1. So these 

2 7 1 ^ 12 

solutions, we see, constitute a fundamental set. 

Next let us assume that not all the first minors of A (A.) vanish when 

o = a 2 = ai ; and let us choose the notation so that the first minor ( 1 J is not 

zero . Then to get a solution associated with a 1} we make the transformation 
y = b 1 x 1 + b 2 x 2 + ' ' ' > y /== ^1 + ^ 2 X ' 2 + ' " ' an d, if possible, deter- 
mine the b's in such a manner that y' = ay. On substituting these equations 
in (1) we get b i [a xx x 1 + a 12 x 2 +•••] + b 2 {a 2X x x r -\- a 22 x 2 + •••]+•••• 

= «[Mi + M! + ' • • ]• 

This equation must hold for all initial values of the x's. Therefore it is an 

identity in them; and the &'s satisfy the following set of equations 

[ («xi — <*)&i + M2+ • • • =0, 
(17) .{ HA + (a 22 — a)l 2 + • • • = 0, 

I • 

In order that these equations have a solution for the &'s not all zero, we get 
as a necessary and sufficient condition, just as in the case of c's in (13), 
A (A) = 0. This condition is satisfied, for we have assumed that when a = %, 
the characteristic equation is satisfied. Furthermore we have chosen the nota- 
tion such that the minor ( 1 ) is not zero when a = a x . Then we get for the 6's 

0j==mj l u j = l,2, • • • 00, 

A) 
where rtij = ^. For convenience we shall take b t = 1. Then the equations 

(1) 
become 

a; 2 =«2i2/i+ («22 — a 21 m 2 )^2 + (a 28 — «2i«i 3 )a ! 3 + - - • , 



(l') : 



fci^/i + ( a <>2 — a ni m 2 )x 2 + (a„ 3 — a ni m z )x s + 
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Next we make the transformation y 2 = d 1 y 1 + d 2 x 2 + • • • , y' = d^ + 
d 2 x' z -\- • • • and determine, if possible, d 2 , d s , • ■ ■ , such that y' = ay 2 + y%- 
Then, as in the foregoing, we get the following, infinite set of equations for 
the determination of the d's 

f (°i — «)^i + « 2l ^2 + « S A +•••=!, 
(17') -> 0+ (a 22 — a 21 m 2 — a)d 2 + ■ • • = 0, 

1 • 

On setting a = — 1/A, we get an equivalent set of equations for the deter- 
mination of the d's, the determinant of which is 



MA) 



1 — j— Actjj A#2x? 

0, 1 + X(a 22 — m 2 a 21 ), 



= (1 + Aa,) 



1 + A(a 22 — a 21 m 2 ),- 
A(a 3 2 — a 21 m 3 ),-- 



By making use of the equations (17), the reader will convince himself that 
A X (A) is identically equal to A (A), for it is obtained from A (A) first by inter- 
changing the rows and columns of A(A) and then multiplying the elements 
of certain rows by certain quantities and adding them to the corresponding 
elements of other rows. Furthermore there is a first minor of A x (A) which is 
not zero when a = %. This minor involves the elements of the first row, 
since all the other first minors vanish. 

The determinant of the coefficients of the equations, omitting the first 
one, is zero for a = %, since when o = a u A(A) has a multiple root A x . There- 
fore we can solve for the ratios of d 2 , d s , • • • . Then let us substitute the 
values of the d's thus obtained in the first equation of (17'). First, we know 
that the sum will converge * ; for the d's are proportional to the first minors 
of A X (A), and the sum of the products of the elements of a row and the cor- 
responding first minors converges. Secondly, since the d's which we have 
determined carry an arbitrary factor, we can finally determine them so that 
the first equation is satisfied. 

Therefore our equations have been reduced to 

(i") y[=^Vi, y' 2 =yi + *iy 2 ,K ==a ?iyi + a sl ) y* + al s 2 s x s + - • • , 

where the a ( { 2) are the transformed at} of the original equations. And on solv- 
ing the systems (1') and (1") and putting in place of y ± and y 2 the cor- 
responding value of x x and x 2 , we see that the solutions of (1) associated with 
a x are 

Xi x = c {l e a i*, Xi 2 = (c i2 + tc il )e a i*„ t= 1, 2, • • • oo. 



F. Eiesz, '-' Les syst&mes d'equations lin&tires a une infinite d'inconnues," p. 34. 
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Then by a precisely similar argument to the foregoing we can show that 

x ix ,Xi 2 ,4>ij, * = 1, 2, • ■ ■ oo, /=3,4, • • • oo, 

constitute a fundamental set. 

When for a given value of a the characteristic equation has a root of higher 

multiplicity, discussions similar to the foregoing must be made. 

We should expect to treat next the case that A(A)= has a root of infinite 
multiplicity. But such a case cannot arise. For we have seen that A (A) is 
an integral transcendental function of A. Therefore it can be expanded in a 
Taylor's series in the neighborhood of any finite A. Suppose that A were a 
finite root of infinite multiplicity. Then A (A) vanishes together with all of 
its derivatives when A = A . Hence A (A) vanishes at every point in the 
neighborhood of A - Then it is identically zero. But this is not true, for 
when A — 0, A (A) = 1. Hence A (A) has no root of infinite multiplicity. 

If a — 0, we can not make the transformation which carries the system 
(12) into the system (13), and the determinant of (IS) diverges. Then we 
do not know whether the equations (12) have a solution for the c§ not all 
zero ; but each special case must be considered as it arises. 

To sum up our results we have 

Theorem IV. If and only if A (A) has 'an infinite number of roots A«, the 
system (1) has a fundamental set of solutions each of whose elements is of 
the form 

x i j = e*j t \l/ij{t), t,; — 1,2, • • • oo, 

where a x = — 1/A» and the tf/ij are polynomials in t of degree at most (n — 1), 
and n is the order of multiplicity of the root Ay. 

Ill 

Fow we shall assume that the 6%j of Theorem I are periodic with the 
period 2ir. We have seen that the system (1) has, as a fundamental set of 
solutions, the set <\>i](t), where 

(18) ^(0) =1, $«,(0) =0, iy^j, »,/, — 1,2, • • • oo. 

Let us make the transformation 

(19) x t = e«*y«, 

where a is an undetermined constant. Then the equations (1) become 

00 

(20) y> + ay i = l On (t)y h i = 1, 2, • • • oo . 
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Since the system 4m constitutes a fundamental set of solutions of (1), any 
solution of (20) can be written 

00 

(21) yi <=e- at %Aj4>ii{t), *=1, 2, • • • oo. 

We now inquire whether it is possible to determine the Aj and a so that 
the yt defined by (21) shall be periodic with the period 2ir. From the form 
of (20) it is evident that necessary and sufficient conditions that the yi be 
periodic with the period 2w are 

(22) yi(2*)— yi(O) = 0, i — 1,2, ■ ■ ■ oo. 
On imposing these conditions on (21), we get 

00 

(23) •%A j [e- 2ajr cj> ij (2w)— <M0)] = 0, i=l,2, • • • oo. 

Then we make the transformation 

(24) e - 2 «T = o-/o- — 1, 

assuming of course that er***^l. Then the equations (23) become 

00 

(25) S'^o-^^tt) + Ai[l +<r(0«(2,r) — 1] = 0, 

i¥*j> * = 1, 2, • ■ • oo. 

In order that these equations have a solution other than that in which the Aj 
are all zero, it is necessary that the determinant of the coefficients be zero. 
On writing the 4m (2tt) simply <j> ih the determinant is 



(26) D(a) — 



l+°"(«ki — 1)> <r<t>i2> 



= 0. 



This determinant, we see from the fact that the determinant of 4> i} converges 
absolutely, is an integral function of a. Equation (26) is called the funda- 
mental equation associated with the period 2tt, and does not admit o- = 
as a root. 

As in the characteristic equation in II, the fundamental equation may 
have no finite root, a finite number of finite roots, or a denumerably infinite 
number of finite roots. 

As in the case of constant coefficients, if the fundamental equation has 
no roots, the system (1) has no solution of the form (19), if the fundamental 
equation has a finite number of finite roots, the system (1) has a finite number 
of solutions of the form (19), and finally if the fundamental equation has an 
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infinite number of finite roots, the system (1) has a fundamental set of solu- 
tions, each of the elements of which is of the form (19). 

The discussion of the form of the solutions follows precisely the same lines 
as that in the case of constant coefficients except when <r is an refold root of 

(26) and not all the (n — fc)th, h > 0, minors of D(a) vanish for <r = o- . 
ISTow let us assume that o- = cr is an n-fold root of (26) and not all the 

first minors of D(<r) vanish. Then there is only one solution of (1) of the 
form x t = e a i*j/ i , where the y% are expressed as in (21). And let us choose 
the notation so that a minor corresponding to the elements of the first column 
is not zero. 

Then we take as a new set of solutions 

(27) x il = e<h t y h> Xij = (j>ij(t), i= 1, 2, • • • oo, j = 2,3, • • • oo. 

This set of solutions can be shown to constitute a fundamental set in a pre- 
cisely similar fashion to that in II. Then we make the transformation 

Xi = e at Zi, i <= 1, 2, • • • oo . 

As above we get 

00 

(28) 2 i = e-«*[4 1 e 0l 'yii + SA^«(f)], »— 1,8, • • • oo. 

1=2 

Necessary and sufficient conditions that the zi be periodic with the period 
2rr are 

Zi(2ir)— 2i(0) =0, t = l, 2, • • • oo. 
On imposing these conditions on (28). we get 

00 

(29) J A 1 [ff*<«-<».'»y« 1 (3») — y H (0)] +SM<M2<r)e- 2 <"r — *«(0)] = 0. 

1=2 

After setting e _2 °' r = cr/<r — 1, the equations (29) become 

00 

(30) A x {l — G/G x )y tx (0) +2'^<7^ + ^[l + o-(^_ 1)]=0, 

1=2 

Tcy^j, k = 2,3, • ■ ■ oo, *=1,2, • • • oo. 



The fundamental equation for the equations (30) is 

(1 — ^1)2/11 (°)> °-<£i2> • • • 
(31) (1 — »/«i)y»(0), 1 +«*(&, — 1), • • • =0. 



Making use of (21) and taking A x = 1, equation (31) becomes 
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00 
(T%Aj<j> 2 j A 2 (<T 1), l+c(022 1), 



:D(<r)— 0. 



So we have 

(32) D(a) = (l-a/a t ) D 1 =(l-a/a 1 ) 



y 1± (0), <r<j) 12 , ■ ■ ■ 

y 21 (0), 1 + <j(4> 22 — 1), ■ ■ -I 



= 0. 



Since o = a 1 is an ra-fold root of D(a) = 0, D t (<T) has the factor (1 — 
a/<r 1 ) n - 1 . Since (27) constitutes a fundamental set, any solution can be 
expressed in the form 

00 

(33) x i = B 1 e a Jyi 1 -\-%Bj<j> 1 j, i=l, 2, • ■ ■ , oo. 

j=2 

Now let us make the transformation, corresponding to (19), to get a 
second solution associated with a x . 

(34) x i2 = e*i*{y i2 -\-ty H ), t" = l,2, • • • oo. 

On imposing the condition that Xi 2 shall satisfy the system (1), we find since 
«° 1 *y*i is a solution, 

00 

( 35 ) y' i2 +<h.yi2 = '2 6i j (t)yj 2 — y ll , i=l, 2, • • • oo. 

From the form of (35) we see that sufficient conditions that the y% 2 shall be 
periodic with the period 2w are 

00 

(36) y i2 (2*)— yt 2 (0) =0 = -2,Bj[e- 2 w<t> ij (2ir) — <M0)] 

— 2^(0) =0. 
On substituting e~ 2a ^ = <t x /<j x — 1, (36) becomes 

00 

(37) — 2^ — 1)^ (0)+ %'B,* i) + BiCl + ff! (*» — 1]=0, Tc^j. 

j=2 



The condition that these equations be consistent is, since o^ =/= 1, 



yu(°)> <*i<t>i2> • • • 

2/2l(0), 1+0-1 (^22 — 1),' 



0. 



In (32) we showed that D x vanishes (n — 1) times when o- = o- 1 . By hypo- 
theses not all the first minors corresponding to the elements of the first column 
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are zero for a^a^. Hence we can solve equations (37) for B 2 , B 3 , • • • , 
terms of 2/11.(0). Consequently in this ease we get a second solution associ- 
ated with a 1} which is of the form (34). 

In a similar manner we can go ahead step by step and get the following 
group of solutions associated with a x 

Xi-L = e a i yi 1} 

x i2 = e<h t ly i2 + ty il ], 



x in = e<h*[y in -f ty in , t + ■ ■ ■ + l/(n — 1) ! t n - 1 y tl '\. 

If <r x is a triple root of D (<x) = 0, such that all first minors are zero, but not 
all the second minors are zero, the solutions associated with % are 

x H = e<h*y H , Xi2 = e^*y i2 , t l3 = e«i'[y 4a + t(y h + y i2 )~]. 

All the sub-cases can be treated, as they arise, by the methods given here. 

IV 

Now let us assume that the system of differential equations (1) satisfy 
the hypotheses of the second existence theorem, and in addition the coefficients 
in the power series expansions of the On are separately periodic with the period 
2ir. That theorem tells us that the solutions of (1) can be written in the form 

00 

(39) out = a *«> (*) ^, i j = 1, 2, • • • 00 . 

And we shall take the initial conditions such that 

00 

(40) *«(0)-S*»(0)^c u , 

whence 

as<« (0) = c th x™ (0) = 0, h = 1, 2, • • • 00. 

where the c tj are constants such that cu = 1, and their determinant is abso- 
lutely convergent. These conditions coupled with the fact that if the deter- 
minant of a set of solutions converges and is not zero when t = 0, it converges 

00 
and is not zero for every value of t for which %6 U converges, show that the 

i=i 
system (40) constitutes a fundamental set of solutions. 

Now we inquire if we can find solutions of (1) of the form 

(41) x i = e^y i> 

where the yi are periodic with the period 2tt, and a is a constant which remains 
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to be determined. After making the transformation (41), the differential 
equations and their solutions become 



(42) 



y > ■ + ayi = s[«„ + 2 6™^]y h i = 1, 8, • • • oo, 

00 00 

y t °-%A,er*W>(t) +Sa?? 1 ) (O^L »— 1,8, • • • oo. 



On imposing the conditions that the y% be periodic with the period %it, viz., 
yi(8ir) — yi (0) = 0, we get 

00 00 

(43) = S ^ [e- 2 «Ta; <0 > (8») — c 4/ + ef*"3te<*] (8*)/^], ♦— 1, 8, • • • oo . 

After setting e- 2ar = <r/o- — 1, the equations (43) become 

00 00 

(44) = a^[ ff {a;<»(2, + a3!«>(a»)^} — a, (a- 1)]=0, 

t = 1, 8, • • • oo . 

To avoid the trivial case where the A) are all zero, we must set the determinant 

00 

(45) D(*,p) = | [>{*<» (2*) +J&>*> (*),»»} — «,(„ — 1)]| = 0, 



which is a condition on the undetermined constant a. Since the x\°) + %x ih) ix h 
are the elements of a fundamental set, and cu = 1, it follows that D(a,p) 
converges absolutely for all finite as. If the fundamental equation is satisfied 
when a = a , it is also satisfied when a = a -f- vV — 1, where v is any 
integer j but all distinct solutions of the differential equations can be obtained 
by taking v = 0, since the ratios of the A } are the same for v = as for v = p. 
And for /t = the equation (45) reduces to 

(46) D .= | «r<« (2») — Cti ( a _ i) | = . 

When /i = 0, we have the case of constant coefficients which we treated 
in II. Now we shall determine the e t] and the a <0) = A" 1 of the solutions as 
we did there. But having the c i} we might also determine the a of the solu- 
tions by means of (46). Since for given initial conditions the solution of the 
differential equations is unique, and since the initial conditions are the same 
in the two cases, the solutions obtained by the former method and the solution 
obtained by means of (46) are the same. 

We recall that the c i} are not all zero if and only if the characteristic 
equation has finite roots. Therefore for every root of the characteristic equa- 
tion there is a corresponding root a- of (46) ; and conversely, for every root 
of (46) there is as corresponding root of the characteristic equation. There- 
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fore if the characteristic equation has no finite roots, then (46) has no finite 
roots, if the former has only a finite number of finite roots, then the latter 
has only a finite number of finite roots, then (46) has an infinite number of 
finite roots which yield independent solutions of the system (1). 

Now since the solution (41) converges for all /*, sufficiently small, includ- 
ing zero, the failure to find an a for /t = is not due to clumsy analytic 
methods, but shows that the system (1) has no solution of the form (41). 
From the theory of implicit functions we know that if (46) is satisfied when 
a = a , then we can solve for a as power series in /t, provided //. is sufficiently 
small. Therefore the existence of a root of the characteristic equation is both 
necessary and sufficient for the existence of a solution of the system (1) of 
the form (41), and the existence of an infinite number of roots of the charac- 
teristic equation is necessary and sufficient for the existence of a fundamental 
set of solutions of (1), each of the elements of which is of the form (41). 

We shall study in detail the case that the characteristic equation has an 
infinite number of roots; and we shall use only those values of a obtained from 
(45) which for ju. = reduce to the values of a obtained from the characteristic 
equation. 

When e~ 2aT = 1, the transformation e' 2aT = <r/<r — 1 can not be made. 
Then as in the case of constant coefficients when a = 0, we can make no gen- 
eral statement about the solution, as the determinant of the A t diverges then. 

V 

Solutions when the a< 0) are distinct and a (0) — a (0) =h mod V — 1. 
The part of (45) which is independent of ju, is 

l_ e aa y (»)T 



D = | aCi^y * — C tj (a — 1) | — | C W | J](l 



„ 2ar 
_ \ ~" ' 1 „-2air 



If (45) were an identity in /*, its roots would be the roots of (46), viz., 
a = o ( y 0) . Let us assume that we have the general case in which it is not an 
identity in ju, and set 



(47) 






a =„(0) +8fc . 






Then we 


get 










(48) 


D („,}»). 


= z> + / iF*(; 


5fc>p.) = 






1 d) 1 (l 


1 — 


_ e- 2 <«<\+y i 


00 

-^nvi-f- 


1 


— <r 2 y 0) ' r 


Mi^ + 1 


— 


e -2(a(») fc+ Sp T 



+ Mft (&,/») j=£k. 
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where Fu {fa, p) is a power series in //. and fa, converging for 

| fa | < oo, | ix. | < P > 0. 

Since by hypothesis no two of the a< 0) differ by an imaginary integer, the 
expansion of (48) as a power series in 8k and fi contains a term in fa of the 
first degree and no term independent of both p. and 8&. Therefore we know by 
the theory of implicit functions that (48) can be solved uniquely for fa as a 
power series of the form 

(49) 8* — /»PfcO»), 

which converges for | p | > but sufficiently small. 

Now we substitute this value of a = af 0) + 8fc in (44), and get an infinite 
number of linear homogeneous equations for the Aj whose determinant con- 
verges and is zero, but the first minors of that determinant are not all zero, 
since by hypothesis, the roots of D = are all distinct and no two differ by 
an imaginary integer. Consequently these equations determine uniquely the 
ratios of the Aj as power series in fi, which converge for p sufficiently small. 
On substituting these ratios in (42) we have the particular solution yuc, 
t = l, 2, • • • oo, expanded as a power series in p. Hence we may write it 

00 

(50) y» = 3 $#>(<)/*'. 

Since the periodicity conditions have been satisfied, 

00 

yilc (t + g,) _ y ik ( t ) - 3[y£ (t + 2*) - y«> (0 V = 0, 
for all ix sufficiently small and for all real t. Therefore 

y% (* + 2 ^ - y% (* ) = o, / = o, i, 2, • • • oo, 

whence it follows that the y 1 **, j = 0, 1, 2, • • • oo, are separately periodic. 

A solution is found in a similar fashion for each a m . 

i 

VI 

Solutions when no two of the a <0) are equal but when a ( 2 0) — a< 0) 

= mod V — 1. 

Suppose that when /t = the characteristic equation has two roots such 
that o ( 2 0) and a^ ' differ by an imaginary integer and that none of the other 
a ( p are congruent to a^ ' mod V — 1. Then we see from (45) 
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(51) 2? (a,/*)- 

i — - — - — i nil — - — -4 — i +^(ft8,) -o, 

where as in (47) we have set a= a (0) + 8 r The term of lowest degree in \ 
alone is found by expanding the first bracket and turns out to be of the second 
degree. To get the terms in /t alone we suppress those involving B t) after 
which we get a factor ju, from each of the first two columns. So we see that 
in general the term of lowest degree in /t alone will be in this case of the 
second degree. Hence we have 

(52) a-KiTi 1 ~ e2Ol(0) ' Tftri 1 -" V) ' 1 

K ' ' n L l — e^Cn^-HO'J 11 L 1 — e 2 <«i <0) + 9 i) ,r J 

+ 8^,(8,, fl )+8 2 F 2 (8 1)ft ) = 0. 

In a similar manner if p of the o (0) are congruent to a< 0) mod V — 1> 
then the term of lowest degree in 8 ± alone is of degree p, and in /t alone it is 
of at least the pth. degree. 

The problem of the form of the solution of (52) is one of implicit func- 
tions. Writing the first terms explicitly we have 

S 2 -J- KuS^ + K odh + terms of higher degree = 0, 

where k u , k 02 , • • • , are constants independent of 8 X and /*. On factoring the 
quadratic terms we get 

(53) (8 t — d x n) (8 X — d 2 n) + terms of higher degree = 0. 

If ^i and d 2 are distinct, there are two solutions, and these have the 
form* 

(54) « u = d lf * + ^P, ( M ) , 8 12 = d& + v?P 2 (/*) , 

where P x and P 2 are power series which converge for n sufficiently small. In 
this case the solutions are found as in V. 

But if d t and d 2 are equal, the character of the solution is in general 
quite different and depends upon the terms of higher degree than the second. 
In general it will be a power series in ±V/»- This case we shall consider 
in detail. 

We see from the form of (53) that the expansion of a x as a power series 
in V7* will contain no term in V /* to the first power, but will have the form 

Suppose that this expansion has been obtained from equation (45). 



* ©hrystal, " Algebra," Vol. 2, pp. 358 ff. 
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Then since a t is not a multiple zero of D, not all the first minors of D are 
zero when o = %. The ratios of the Aj will be determined from (44). If 
juA is a non-vanishing first minor corresponding to an element in the first 
column of J), it follows from the form of (45), remembering that we have 
x<p(t) = Cije a( ° u , that solving (44), we get 

A^^AvA^^A^ ; = 3,4, ..-, oo, 
fiA /xA 

where A = A (0) + A'^y/* _|_ a' 1 '/* + , 

A,- = A <0) + A<V2y/2 _j_ A » )/t _|_ 

On substituting these series for the Aj in (42) we find that the yi ± are devel- 
opable as series of the form 

yn = C + vT^' 2 + n> + » - 1, a, ■ • • oo. 

So we see that in general the y« x carry terms in Vp" although the term in V/* 
is absent in the expansion for %. 

However, if all the first minors corresponding to the elements of the first 
column are zero, and if there is a first minor distinct from zero corresponding 
to the elements of the second column, the results are precisely the same. But 
suppose that all the first minors corresponding to the elements of both the 
first and second columns are zero. Then suppose that a first minor corres- 
ponding to an element of the fcth column is not zero. Then it follows from 
the form of (45) that when o = o;, it will carry the factor p 2 ; and let this 
minor be denoted by n 2 A. Then solving (44) we get 



A r"l A A ***** A A ** *■*' A o A 

^i^ — ij, A s = -~ -A k , A j = - TZ -Aj t , ; = 3,4,- • • oo. 



^ A A _ ^ A A ^ 

/rA ti 2 A /t 2 A 

where A 15 A 2 , A,-, do not in general vanish when ai = 0. It follows from the 
first two equations that Ak must carry ftasa factor, since the Aj are finite 
for ii = 0. Hence in this case the y< t have the same form as before. Simi- 
larly the yi 2 have the same properties. 

The solutions associated with a< 0) , a< 0) , • • • , are found as in the pre- 
ceding case. If there are several groups of a (0) in which these congruences 
exist the discussion must be made for each one separately. 

VII 

Solutions when o <0) is a multiple root. 

Now suppose that two of the a< 0) are equal and only two, and that there 
are none of the congruences treated in VI. Let us choose the notation so 
that a'"' = a ( 2 0) . Prom our work in the theory of linear differential equations 
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in infinitely many variables in constant coefficients and from IV, the solutions 
are of the form 

(54) a;i=S4 J -[c^e°/° ) * + S i a;«y] or 

(55) Xi = 4 1 [e ll «*< , >' + a as«y] + A 2 {{a 2 + fc«i)e*<*>* + 3s<*V*] 

00 00 

After setting e' 2av = a/a — 1, the fundamental equation becomes either 

00 

(56) D(a, j») = I [a(c tie 2 «/°>' +JU*> (fcr),.*) - e ti (a— 1)] | = or 

00 

(57) 0(a, M ) = I [»{^^ + S*<» (8») M *} _ Cil (a-l)], 



[*{(«*» + 27rc <x )e 2 M ^ + s J.GB (a,) M *} _ c« s (a — 1)], 



= 0. 



For /x=0 both of these, by the theory of infinite determinants, reduce to 
(58) D — | CJi |[l — a(l — e*^"")]» n [1 — „(l — e 2«/»)T)]. 



3=3 



As before we set a = a^ + 8 W and expand as a power series in 8 X . Then we 
see that the term of lowest degree in S t alone is of the second. When the 
determinant is of the form (56) with a<°> = a<°\ the term of lowest degree 
in n alone is of the second in general. Then we have the same form as in VI. 
But if the determinant is of the form (57), the term in n alone is in general 
of the first degree. In the former case we have a consideration similar to 
that in VI; in the latter case, in general the solutions for 8 X are of the form 
(59) 8ii = Al i/ 2j p (/tl /2) 

where P is a power series in y/J and contains a term independent of p. The 
discussion of the special eases is made just as in VI. On substituting these 
expansions for a = a< i '»+8 l in (44) we solve for the A, as power series in 
\V These Aj substituted in (42) give y H and y H as power series in V£ 
If p of the a<°> are equal, then for these roots the expansions of D starts 
with 8p as the term of lowest degree in S ± alone, and except in the special cases 
corresponding to those mentioned in the foregoing, the term in p. alone is of 
the first degree. Consequently in general for a (0 > = a m ■ ■ ■ = o (0) we have 

12 p > 

**i-WP(W), y— o,i, • • -, p_ i, 

where « is a pth root of unity. 
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VIII 

Solutions when there are equalities and congruences among the a (0) . 

Suppose that two of the a (0> , for example a^ and o^' are equal, and that 
a third one, say a ( 3 0) , differs from a* * by an imaginary integer. Furthermore 
we shall assume that there are no other equalities or congruences among the 
O co)_ ij wo cases ar j se here: (a) the solutions are of the form (54) with 
o < 0) = a^ 0) ; (b) the solutions are of the form (55) 

Case (a). In this case we have 

D =.\c t ,\ [1 — <t(1 — e 2a i (0 >"-)] 3 n [i_ ff (i_ e «ajc>T)] # 

In setting a = a'"' + S lf we find that the term in D of lowest degree in S t 
alone is of the third. To get the term in D of lowest degree in /x alone we set 
S x = 0. Then it becomes evident at once that each of the first columns carry 
/iasa factor, while the remaining ones do not. Consequently the term of 
lowest degree in /* is of the third degree at least. Furthermore since the first 
three columns vanish when //. = S ± = there are no terms of lower degree 
than the third in B ± and /t. Hence in general we see D of the form 

(60) D = 8» + y 21 S^ + y 03 ^ + • • • = 0. 

The problem is now one of implicit functions. The details of the special oases 
must be treated as they arise. However, we make the general statement that 
since the roots of the cubic terms of (60) set equal to zero are in general 
distinct, it follows from the theory of implicit functions that the three values 
of 8 X are in general expansible in integral powers of p. 
Case (b). In this case we have 

D =[l — tr(l — e 2 ^*)] 3 n [1 — <r(l — e*»/Or)— (J. 

On introducing 8 X as before, we find that the term of lowest degree in 8 X alone 
is of the third degree. But when the terms involving /x are retained in D, only 
the first and third columns vanish when /t = 8 1 = 0, and consequently the 
expansion of D will contain a term in ju. 2 alone. Furthermore, since the first 
and third columns vanish for /a = S t = 0, there will be no terms of degree 
lower than the second in /* and 8^ Hence in general D has the form 

(61) D = 8^ + Vll 8^ + y 02 /* 2 + • • • = 0. 

In the general case in which y n and y 02 are not zero, there is one solution in 
integral powers of p. and two in powers of V7'- 

When the roots o< 0) have higher multiplicities and more congruences 
among them, we make a similar discussion. 



